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ABSTRACT 

The influence of mass inliomogeneities on tlie angular size-redshift test 
is investigated for a large class of flat cosmological models driven by dark 
energy plus a cold dark matter component. The results are presented in 
two steps. First, the mass inhomogeneities are modeled by a generalized 
Zeldovich-Kantowski-Dyer-Roeder (ZKDR) distance which is characterized by 
a smoothness parameter a{z) and a power index 7, and, second, we provide 
a statistical analysis to angular size data for a large sample of milliarcsecond 
compact radio sources. By marginalizing over the characteristic angular size 
I, fixing = 0.26 and assuming a Gaussian prior on cu, i.e., cu — —1 ± 0.3, 
the best fit values are cu — —1.03 and a — 0.9. By assuming a Gaussian prior 
on the matter density parameter, i.e., Qm = 0.3 ± 0.1, the best fit model for a 
phantom cosmology with lu = —1.2 occurs at Qm = 0.29 and a = 0.9 when we 
marginalize over the characteristic size of the compact radio sources. The results 
discussed here suggest that the ZKDR distance can give important corrections 
to the so-called background tests of dark energy. 

Subject headings: cosmology: observations - ZKDR distance - phantom dark 
energy 
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1. Introduction 

An impressive convergence of recent astronomical observations are suggesting that 
our world behaves hke a spatially flat scenario dominated by cold dark matter (CDM) plus 
an exotic component endowed with large negative pressure, usually named dark energy 
(Perlmutter et al. 1998; Efstathiou et al. 2002; Riess et al. 1999, 2004; Allen et al. 2004; 
Astier et al. 2006; Riess et al. 2007). In the framework of general relativity, besides the 
cosmological constant, there are several candidates for dark energy, among them: a vacuum 
decaying energy density, or a time varying A{t) (Ozer & Taha 1986, 1987; Bertolami 1986; 
Preese et al. 1987; Carvalho et al. 1992; Lima & Maia 1994; Lima & Trodden 1996; Lima 
1996; Torres & Waga 1996; Overduin & Cooperstock 1998; Cunha & Santos 2004; Shapiro 
et al. 2005, Costa et al. 2007), the so-called "X-matter" (Turner & White 1997; Chiba et 
al. 1997, Alcaniz & Lima 1999, 2001, Cunha et al. 2003, D^browski 2007), a relic scalar 
field (Peebles & Ratra 1998; Caldwell et al. 1998; Ulam et al. 2004), and a Chaplygin 
Gas (Kameshchik et al. 2001; Bihc et al. 2002; Bento et al. 2002; Alcaniz & Lima 2005). 
Some recent review articles discussing the history, interpretations, as well as, the major 
difficulties of such candidates have also been published in the last few years (Padmanabhan 
2003; Peebles & Ratra 2003; Lima 2004; Turner & Huterer 2007). 

In the case of X-matter, for instance, the dark energy component is simply described by 
an equation of state = uip^. The case ui — —1 reduces to the cosmological constant, and 
together the CDM defines the scenario usually referred to as "cosmic concordance model" 
(ACDM). The imposition a; > — 1 is physically motivated by the classical fiuid description 
(Hawking & Ellis 1973). However, as discussed by several authors, such an imposition 
introduces a strong bias in the parameter determination from observational data. In order 
to take into account this difficulty, superquintessence or phantom dark energy cosmologies 
have been recently considered where such a condition is relaxed (Faraoni 2002; Caldwell 
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et al. 2003; Gonzales-Diaz 2003, Santos & Alcaniz 2005, Linder 2007). In contrast to 
the usual quintessence model, a decoupled phantom component presents an anomalous 
evolutionary behavior. For instance, the existence of future curvature singularities, a 
growth of the energy density with the expansion, or even the possibility of a rip-off of 
the structure of matter at all scales are theoretically expected (see, however, Alcaniz & 
Lima 2004; Gonzalez-Diaz & Siguenza (2004), de Freitas Pacheco & Hovarth (2007) for 
a thermodynamic discussion). Although possessing such strange features, the phantom 
behavior is theoretically allowed by some kinetically scalar field driven cosmology (Chiba et 
al. 2000), as well as, by brane world models (Shani & Shtanov 2002, 2003, Wu et al. 2007), 
and, perhaps, more important to the present work, a PhantomCDM cosmology provides 
a better fit to type la Supernovae observations than does the ACDM model(Alam et al. 
2003; Choudury & Padmanabhan 2004; Astier et al. 2006). Many others observational and 
theoretical properties of phantom driven cosmologies have also been successfully confronted 
to standard results (see, for instance, Alcaniz 2004; Piao & Zhang 2004; Choudury and 
Padmanabhan 2005; Perivolaropoulos 2005; Nesseris & Perivolaropoulos 2007). 

In this context, one of the most important tasks for cosmologists nowadays is to 
confront different cosmological scenarios driven by cold dark matter (CDM) plus a given 
dark energy candidate with the available observational data. As widely known, a key 
quantity for some cosmological tests is the angular distance- redshift relation, Da{z), which 
for a homogeneous and isotropic background, can readily be derived by using the Einstein 
field equations for the Priedmann- Robertson- Walker (FRW) geometry. Prom Da{z) one 
obtains the expression for the angular diameter 6{z) (see section 3) which can be compared 
with the available data for different samples of astronomical objects (Gurvits et al. 1999; 
Lima & Alcaniz 2000, 2002; Gurvits 2004; Alcaniz & Lima 2005). 

Nevertheless, the real Universe is not perfectly homogeneous, with light beams 
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experiencing mass inhomogeneities along their way. Actually, from small to intermediate 
scales (< lOOMpc), there is a lot of structure in form of voids, clumps and clusters which 
is probed by the propagating light. Since the perturbed metric is unknown, an interesting 
possibility to account for such an effect is to introduce the smoothness parameter a which 
is a phenomenological representation of the magnification effects experienced by the light 
beam. From general grounds, one expects a redshift dependence of a since the degree 
of smoothness for the pressureless matter is supposed to be a time varying quantity 
(Linder 1988). When a — 1 (filled beam), the homogeneous FRW case is fully recovered; 
a < 1 stands for a defocusing effect while a = represents a totally clumped universe 
(empty beam). The distance relation that takes these mass inhomogeneities into account 
was discussed by Zeldovich (1964) followed by Kantowski (1969), although a clear-cut 
application for cosmology was given only in 1972 by Dyer & Roeder (many references may 
be found in the textbook by Schneider, Ehlers & Falco 1992; Kantowski 2003). Many studies 
involving the ZKDR distances in dark energy models have been published in the literature. 
Analytical expressions for a general background in the empty beam approximation {a — 0) 
were derived by Sereno et al. (2001). By assuming that both dominant components may 
be clustered they also discussed how the critical redhift, i.e., the value of z for which 
Da{z) is a maximum (or Q{z) minimum), and compared to the homogeneous background 
results as given by Lima & Alcaniz (2000), and, further discussed by Lewis & Ibata (2002), 
and Araiijo & Stoeger (2007). More recently, Demianski et al. (2003), derived an useful 
analytical approximate solution for a clumped concordance model (ACDM) valid on the 
interval < 2; < 10. Additional studies on this subject involving time delay (Giovi & 
Amendola 2001; Lewis & Ibata 2002) gravitational lensing (Kochanek 2002; Kochanek & 
Schechter 2003) or even accelerated models driven by particle creation (Campos & Souza 
2004) have also been considered. 

Although carefully investigated in many of their theoretical and observational aspects. 
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an overview in the literature shows that a quantitative analysis on the influence of dark 
energy in connection with inhomogeneities present in the observed universe still remains to 
be studied. Recently, the ZKDR distance was applied for the 9{z) statistics with basis on 
a ACDM cosmology with constant a ( Alcaniz et al. 2004) . It was concluded that the best 
fit model occurs at Qm — 0.2 and a — 0.8 whether the characteristic angular size I of the 
compact radio sources is marginalized. 

In this paper, we focus our attention on X-matter cosmologies with special emphasis 
to phantom models (a; < — 1) by taking into account the presence of a clustered cold dark 
matter. The mass inhomogeneities will be described by the ZKDR distance characterized 
by a smoothness parameter a{z) which depends on a positive power index 7. The 
main objective is to provide a statistical analysis to angular size data from a large 
sample of milliarcsecond compact radio sources distributed over a wide range of redshifts 
(0.011 < z < 4.72) whose distance is defined by the ZKDR equation. As an extra bonus, 
it will be shown that a pure CDM model (Qm = 1) is not compatible with these data even 
for the empty beam approximation (a = 0). The manuscript is organized as follows. In 
section 2 we derive the ZKDR equation. We also provide some arguments for a locally 
nonhomogeneous Universe where the homogeneous contribution of the dark matter obeys 
the relation = apo^Pml PoV where 7 is a positive number. In section 3 we analyze the 
constraints on the free parameters a, 7 and VLm- We end the paper by summarizing the 
main results in section 4. 



2. The Extended ZKDR Equation 



Let us now consider a flat FRW geometry (c = 1) 




(1) 
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where R{t) is the scale factor. Such a spacetime is supported by the pressureless CDM fluid 
plus a X-matter component of densities pm and p^, respectively. Hence, the total energy 
momentum tensor, T^^ = T^^(^m) + T^^(^x), can be written as 

T'^" - [pM + (1 + uj)p,]UW - ajp^g''", (2) 

where Lf^ = 6^ is the hydrodynamics 4- velocity of the comoving volume elements. In this 
framework, the Einstein Field Equations (EFE) 

G^" = Rf"" - -g^^R = SttCT'^^ (3) 

take the following form: 

/? _ r n n . .1 

(4) 



R ^ -2^" 



(5) 



R' ' ' R 

where an overdot denotes derivative with respect to time and Hg = 100/iKms~^Mpc~^ is 
the Hubble parameter. By the flat condition, fl^ — i — f^M, is the present day dark energy 
density parameter. As one may check from (2)-(5), the case u; — —1 describes effectively 
the favored "cosmic concordance model" (ACDM). 

On the other hand, in the framework of a comformally flat FRW metric, the optical 
scalar equation in the geometric optics approximation reads (Optical shear neglected) 



Va" + ^R^^ki^k^yfA = 0, (6) 

where A is the beam cross sectional area, plicas means derivative with respect to the 
affine parameter describing the null geodesies, and A;^ is a 4-vector tangent to the photon 
trajectory whose divergence determines the optical scalar expansion (Linder 1988; Giovi 
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& Amendola 2001; Demianski et al. 2000, Sereno et al. 2001). The circular frequency of 
the hght ray as seen by the observer with 4- velocity C/" is a; = C/°fca, while the angular 
diameter distance, Da, is proportional to \/A (Schneider, Ehlers & Falco 1992). 

As widely known, there is no an acceptable averaging procedure for smoothing out local 
inhomogeneities. After Dyer & Roeder (1972), it is usual to introduce a phenomenological 
parameter, a{z) — 1 — ^ , called the "smoothness" parameter. For each value of z, 
such a parameter quantifies the portion of matter in clumps (pd) relative to the amount of 
background matter which is uniformly distributed (pm)- As a matter of fact, such authors 
examined only the case for constant a, however, the basic consequence of the structure 
formation process is that it must be a function of the redshift. Combining equations (2), 
(3) and (6), after a straightforward but lengthy algebra one finds that the angular diameter 
distance, Da{z), obeys the following differential equation 

(l + -)^^^+(l + -)6^ + W. = 0, (7) 

which satisfies the boundary conditions: 

Da (0) = 0, 



(8) 



dDA I 



1. 



The functions J-', Q and Ti in equation (7) read 



Q = \^M{^^zf^'^^^{\-^M){^^zf^-^'^ 



3(a; + 1 



1-0m)(1 + ^)'("+'^. (9) 
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Fig. 1. — The smoothness parameter as a function of the redshift for some selected values 
of Po and 7. All curves approach the filled beam result (a = 1) at high redshifts regardless 
of the values of Po and 7. Note that Po determines cto = a{z = 0). For a given Po the 
curves starts at the same point but the rate approaching unit (filled beam) depends on the 
7 parameter. 

The smoothness parameter a{z), appearing in the expression of Ti, assumes the form below 
(see Appendix A for a detailed discussion) 

Po{l + zY^ 

<^(-2) = TTT. (10) 

where Po and 7 are constants. Note that the fraction ao = /3o/(l + Po) is the present day 
value of a{z). In Fig. 1 we show the general behavior of a{z) for some selected values of Po 
and 7. 

At this point, it is interesting to compare Eq. (7) together the subsidiary definitions 
(8)-(10) with other treatments appearing in the literature. For 7 = (constant a) and 
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o; = -1 (ACDM) it reduces to Eq. (2) as given by Alcaniz et al. (2004). In fact, for a; = -1 
the function V, is given by 7^ = ^^2^(1 + z^. A more general expression for ACDM model 
(by including the curvature term) has been derived by Demianski et al. (2004). As one may 
check, by identifying u = m/3 — 1, our Eq. (7) is exactly Eq.(lO) presented by Giovi & 
Amendola (2001) in their time delay studies (see also Eq. (2) of Sereno et al. (2002)). It 
is worth notice that in this paper the a parameter is greater than unity. This means that 
the light rays are demagnified along the path (for a < 1 the light rays are magnified). In 
addition, the a parameter may also depend on the direction along the line of sight (for a 
discussion of such effects see Wang 1999, Sereno et al. 2002). 

Let us now discuss the integration of the ZKDR equation with emphasis in the so-called 
phantom dark energy model (a; < —1). In what follows, assuming that a; is a constant, we 
have applied for all graphics a simple Runge-Kutta scheme (see, for instance, the rksuite 
package from www.netlib.org). 

In Figure 2 one can see how the equation of state parameter, uj, affects the angular 
diameter distance. For fixed values of Qm — 0.3, Po — 0.5 and 7 = 0, all the distances 
increase with the redshift when cu diminishes and enters in the phantom regime {cu < —1). 
For comparison we have also plotted the case for ACDM cosmology (a; = —1). 

In Fig. 3 we show the effect of the 7 parameter on the angular diameter distance 
for a specific phantom cosmology with u = —1.3, as requested by some recent analyzes 
of Supernovae data (Riess 2004, Perivolopoulus 2004). For this plot we have considered 
Po — 0.5. As shown in Appendix A, /3o = {Ph/ Pci)z=o, is the present ratio between the 
homogeneous (ph) and the clumped (pd) fractions. It was fixed in such a way that ao 
assumes the value 0.33. Until redshifts of the order of 2, the distance grows for smaller 
values of 7, and after that, it decreases following nearly the same behavior. 

In Fig. 4 we display the influence of the Po parameter on the angular diameter distance 
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Fig. 2. — Angular diameter distance for a flat FRW phantom cosmology. The curves display 
the effect of the equation of state parameter for f3o = 0.5 and 7 = 0. The thick curve 
corresponds to the ACDM model. Note that for a given redshift, the distances always 
increase for u beyond the phantom divide line {u < —1). 

for two distinct sets of 7 values. The cosmological framework is defined = 0.3 and the 
same equation of state parameter u = —1.3 (phantom cosmology). For each branch (a 
subset of 3 curves with fixed 7) the distance increases for smaller values of as should be 
expected. 

3. ZKDR distance and Angular Size Statistics 

As we have seen, in order to apply the angular diameter distance to a more realistic 
description of the universe it is necessary to take into account local inhomogeneities in the 
distribution of matter. Similarly, such a statement remains true for any cosmological test 
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Fig. 3. — Effects of the 7 parameter on the angular diameter distance. For all curves we 
fixed oj = —1.3, Po = 0.5 and Qm = 0.3. Note that the distances increase for smaller values 
of 7. 

involving angular diameter distances, as for instance, measurements of angular size, 0{z), 
of distant objects. Thus, instead of the standard FRW homogeneous diameter distance one 
must consider the solutions of the ZKDR equation. 

Here we are concerned with angular diameters of light sources described as rigid rods 
and not isophotal diameters. In the FRW metric, the angular size of a light source of proper 
length / (assumed free of evolutionary effects) and located at redshift z can be written as 



Da(z] 



(11) 



where i = 100/h is the angular size scale expressed in milliarcsecond (mas) while / is 
measured in parsecs for compact radio sources (see below). 



Let us now discuss the constraints from angular size measurements of high z objects 



- 13 - 



< 



0.0 



0.4 - 



0.3 - 



0.2 - 



0.1 








2 



4 6 



8 



10 



Redshift (z) 



Fig. 4. — Influence of tlie /3o parameter on tlie angular diameter distance for ^Im = 0.3 and 
uj = —1.3. The curves are separated in two sets corresponding to the values of 7 = 0.5,0.9 
as indicated in the box. As expected, both sets present the same behavior at low redshifts. 

on the cosmological parameters. The present analysis is based on the angular size data 
for milliarcsecond compact radio sources compiled by Gurvits et al. (1999). This sample 
is composed by 145 sources at low and high redshifts (0.011 < z < 4.72) distributed into 
12 bins with 12-13 sources per bin (for more details see Gurvits et al. 1999). In Figure 5 
we show the binned data of the median angular size plotted as a function of redshift z to 
the case with 7 = and some selected values of VLm and = /3o/(l — /?o) = constant. As 
can be seen there, for a given value of VLm the corresponding curve is slightly modified for 
different values of the smoothness parameter a. 

Now, in order to constrain the cosmic parameters, we first fix the central value of the 
Hubble parameter obtained by the HST key project ifo = 72 ± 8 km.s^^.Mpc^^ (Freedman 
et al. 2001). Note that this value is greater that the recent determination by Sandage and 
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Fig. 5. — Angular size versus redshift according to the ZKDR distance. Curves for Qm = 0.3, 
7 = and different values of uj are shown. The data points correspond to 145 compact radio 
sources binned into 12 bins (Gurvits et al. 1999). For comparison the filled beam ACDM 
has been included. 

collaborators (see astro-ph/0603647), and it is in accordance with the 3 years release of the 
WMAP team. Following standard lines, the confidence regions are constructed through a 
minimization 



where 9{zi, I, u, a) is defined from Eq. (7) and 9oi are the observed values of the angular 
size with errors cTj of the ith bin in the sample. The confidence regions are defined by the 
conventional two-parameters levels. In this analysis, the intrinsic length /, is considered 
a kind of "nuisance" parameter, and, as such, we have also marginalized over it. 

In Fig. 6 we show confidence regions in the u — a plane fixing Qm = 0.263, and 
assuming a Gaussian prior on the u parameter, i.e., u = —1 ± 0.3 (in order to accelerate 




(12) 
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Fig. 6. — Confidence regions in the uj — a plane according to the sample of angular size data 
by Gurvits et al. (1999) and fixed = 0.263 as shown in panel. The confidence levels of 
the contours are indicated. The point "x" marks the best fit values, u = —1.03 and a = 0.90. 

the universe). The "x" indicates the best fit model that occurs at a; = —1.03 and a ~ 0.9. 

In Fig. 7 the confidence regions are shown in the Qm ~ <^ plane. We have now assumed 
a Gaussian prior on VLm-, i-e., = 0.3 ±0.1 from the large scale structure. From Figs. 6 
and 7, it is also perceptible that while the parameters uj and VLm are strongly restricted, 
the entire interval of a is still allowed. This shows the impossibility of tightly constraining 
the smoothness parameter a with the current angular size data. This result is in good 
agreement with the one found by Lima & Alcaniz (2002) where the same data set were used 
to investigate constraints on quintessence scenarios in homogeneous background, and is also 
in line with the one obtained by Barber et al. (2000) who argued in favor of Uo = a{z = 0) 
near unity (see also Alcaniz, Lima & Silva 2004 for constraints on the ACDM model). 
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Fig. 7. — Confidence regions in the Qm — a plane according to the sample of angular size 
data by Gurvits et al. (1999). For a phantom cosmology with u = —1.2, the confidence 
levels of the contours are indicated. As in Fig. 6, the "x" also points to the best fit values 
shown in the panel. 

4. Summary and Concluding Remarks 



All cosmological distances must be notably modified whether the space-time is filled 
by a smooth dark energy component with negative pressure plus a clustered dark matter. 
Here we have addressed the question of how the angular diameter distance of extragalactic 
objects are modified by assuming a slightly inhomogeneous universe. The present treatment 
complements our previous studies (Alcaniz & Lima 2000, 2002) by considering that the 
inhomogeneities can be described by the Zeldovich-Kantowski-Dyer-Roeder distance (in 
this connection see also, Giovi & Amendola 2001; Lewis & Ibata 2002; Sereno et al. 2001; 
Demianski et al. 2003). The dark energy component was described by the equation of state 
= ujp^. A special emphasis was given to the case of phantom cosmology {uj < —1) when 
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the dominant energy condition is violated. The effects of the local clustered distribution of 
dark matter have been described by the "smoothness" phenomenological parameter a{z), 
and a simple argument for its functional redshift dependence was given in the Appendix A 
(see also Figure 1). 

The influence of the dark energy component was quantified by considering the angular 
diameters for sample of milliarcsecond radio sources (Fig. 5) as described by Gurvits et al. 
(1999). By marginalizing over the characteristic angular size I and assuming a Gaussian 
prior on the matter density parameter, i.e., Qm = 0.263 ± 0.07(stat) ±0.032 (sys), the best 
fit model occurs at a; = —1.03 and a — 0.9. This phantom model coincides with the central 
value recently determined by the Supernova Legacy Survey (Astier et. al. 2006). On the 
other hand, fixing co = —1.3 and assuming a Gaussian prior for D,m, that is, — 0.3 ±0.1, 
we obtained the best fit values {Qm = 0.29, a = 0.9). 

Finally, it should be stressed that measurements from the angular size combined with 
the ZKDR approach may provide an important and more rigorous cosmological test in the 
near future. However, it is necessary a statistical study for determining the intrinsic length 
of the compact radio sources in order to improve the present results. 



A. On the redshift dependence of a{z) 

In this Appendix we discuss the functional redshift dependence of the smoothness 
parameter, a{z), adopted in this work. By definition 

= 1 _ (Al) 

where pd denotes the clumped fraction of the total matter density, pm, present in the 
considered FRW type Universe. This means that the ratio between the homogeneous 
(ph) and the clumped fraction can be written as Ph/pd — ck{z)/1 — a{z). How this ratio 
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depends on the redshift? In this concern, we first remember that a{z) hes on the interval 
[0,1]. Secondly, in virtue of the structure formation process, one expects that the degree of 
homogeneity must increase for higher redshifts, or equivalently, the clumped fraction should 
be asymptotically vanishing at early times, say, for z > 100. This means that a{z) — > 1 
at high z. At present, {z — 0), this fraction may have an intermediate value, say, /3o- In 
addition, it is also natural to suppose that the redshift dependence of the total matter 
density, pm, must play an important role in the evolution of their fractions. In this way, for 
the sake of generality, we will assume a power law 

p, ^ a{z) ^ p 
Pa 1 - a{z) po 

where /3o = {Ph/Pci)z=o ^-nd 7 are dimensionless numbers. Finally, inserting Pm{z), and 
solving for a{z) we obtain: 

which is the expression adopted in this work (see Eq. (9)). 

As one may check, for positive values of 7, the smoothness function (A. 3) has all the 
physically desirable properties above discussed. In particular, the limit for high values of z 
does not depend on the values of jSo and 7 (both of the order of unity). Note also that if 
the clumped and homogeneous portions are contributing equally at present {Po = 1), we see 
that a{z — 0) — 1/2 regardless of the value of 7. Figure 1 display the general behavior of 
a{z) with the redshift for different choices of Po and 7. The above functional dependence 
should be compared with the other ones discussed in the literature (Linder 1988, 1998; 
Campos & de Souza 2004 and Refs. therein). One of the most interesting features of (A. 3) 
is that its validity is not restricted to a given redshift interval. 
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